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A general approach for constructing chaotic synchronized dynamical systems is discussed that is 


based on a decomposition of given systems into active and passive parts. 


As a possible application 


we consider an improved encoding method where the information signal is injected into the dynamical 
system of the transmitter. Furthermore, we show how to design in a systematic way high-dimensional 
synchronized systems that may be used for efficient hyperchaotic encoding of information. 


PACS numbers: 05.45.+b, 43.72.+q, 47.52.+j 


Synchronization of periodic signals is a well-known 
phenomenon in physics, engineering, and many other sci- 
entific disciplines. Recently, synchronization of chaos [1] 
has aroused much interest in light of its potential appli- 
cations. In particular, the use of chaotic synchronization 
in communication systems has been investigated by sev- 
eral authors [2-11]. There, an information signal con- 
taining a message is transmitted using a chaotic signal as 
a broadband carrier, and the synchronization is necessary 
to recover the information at the receiver. Different im- 
plementations of this basic idea have been suggested. For 
example, in Refs. [2—5] the information signal is added to 
the chaotic signal and in Refs. [2,6] a parametric modula- 
tion is used for the transmission of digital signals. Other 
approaches to use chaos for the purpose of communication 
include controlling techniques to encode binary informa- 
tion [12] and methods that make use of the quick decay of 
the correlation function for chaotic signals [13]. 

In this Letter we discuss a new approach [8—11] for 
constructing (chaotic) synchronized systems that may be 
viewed as a generalization of the method introduced by 
Pecora and Carroll [1]. This approach and two exam- 
ples for illustration are presented in the first part of the 
Letter. In the second part we apply a new method for 
encoding messages using chaotic dynamics [8-11]. In 
contrast to most of the schemes proposed in the literature 
until now we consider cases where the information signal 
drives the dynamical system that is used in the transmitter. 
The scalar signal which is transmitted from the transmit- 
ter to the receiver is a function of the transmitter state 
variables and the information signal. If the receiver syn- 
chronizes with the transmitter, the information signal can 
be recovered exactly, i.e., without the reconstruction error 
that typically occurs with other encoding methods based 
on synchronization [2,3]. Furthermore, it turns out that 
this modulation technique not only yields a transmission 
without errors but also a more secure encoding. Finally, 
in the third part of the Letter we show that the new syn- 
chronization method can be used to construct systemati- 
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cally high-dimensional synchronized systems using low- 
dimensional systems as building blocks. This possibility 
is, for example, very useful for the design of communica- 
tion systems that are based on hyperchaotic signals. 

The new synchronization method is based on the fact 
that it is possible to consider more general decompositions 
of a given dynamical system 


z = F(z) (1) 


than the decomposition into subsystems proposed by 
Pecora and Carroll [1]. Starting from a chaotic autono- 
mous system, for example, one can always formally 
rewrite it in different ways as a nonautonomous system 


x = f(x, s(t) (2) 
with some driving s(t) = h(x) or s = A(x,s). Let 
y = f(y, s(t) (3) 


be a copy of the nonautonomous system that is driven by 
the same signal s(t). If the differential equation for the 
difference e = x — y, 


èe = f(x,s) — f(y,s) = f(x,s) — f(x — e,s), (4) 


possesses a stable fixed point at e = 0 then there exists for 
the systems (2) and (3) a synchronized state x = y that is 
stable. This can be proved using stability analysis of the 
linearized system for small e or using (global) Lyapunov 
functions. In general, however, the stability has to be 
checked numerically using the fact that synchronization 
occurs if all conditional Lyapunov exponents [1] of the 
nonautonomous system (2) are negative. In this case 
system (2) is a passive system that tends to a fixed point 
when not driven. Therefore, we call the decomposition 
given by h and f an active-passive decomposition (APD) 
of the original dynamical system (1). The fact that all 
conditional Lyapunov exponents of (2) are negative does 
not exclude chaotic solutions. 


© 1995 The American Physical Society 


VOLUME 74, NUMBER 25 


PHYSICAL REVIEW LETTERS 


19 JUNE 1995 


To illustrate this synchronization scheme we consider 
different active-passive decompositions of the well-known 
Lorenz model. In the first example we choose 


x; = —10x, + s(t), 
X2 = 28x; — X2 — X1X3, (5) 
X3 = x1X2 — 2.666x3, 
with 
s(t) = h(x) = 10x2. 


To estimate the temporal evolution of the difference 
e =x — y of the states of the two systems (2) and 
(3) we note first that the difference e} = x, — yı of 
the first components converges to zero, because é; = 
—10e,. Therefore, the remaining two-dimensional system 
describing the evolution of the differences e2 = x. — y2 
and e3 = x3 — y3 can for the limit t — © be written as 


e2 = — l2 ~ X123, 
e3 = X162 — 2.666e3 . 


Using the Lyapunov function L = eż + eå one can show 
that L = —2(e5 + 2.6663) < 0. This means, that the 
synchronization is globally stable and occurs for all 
types of driving signals s(t). The conditional Lyapunov 
exponents of this decomposition are given by A; = 
—1.805, Ax = —1.861, and A; = —10 with respect to the 
natural logarithm. 

Our second example is essentially a decomposition of 
the Lorenz system into subsystems as suggested by Pecora 
and Carroll [1]: 


II 


xı = 28s(t) — xı — s(t)x2, 
X2 = s(t)xy — 2.666x> , 


where the function s(t) is now given by the additional 
differential equation 


s = A(x,s) = 10(x; — s). 


Using similar arguments as for the first example (5) it is 
easy to see that the difference e = x — y converges to 0 
for all s(t) [14]. The conditional Lyapunov exponents are 
A; = —1.796 and Az = —1.870. 

This example shows that the Pecora-Carroll method for 
constructing synchronized systems is included in the more 
general APD approach. However, our numerical simula- 
tions [11] indicate that using the APD any typical chaotic 
system can be used to implement synchronization using 
a one-dimensional drive in many ways. The freedom to 
choose the function h that defines the driving signal there- 
fore leads to a large flexibility in applications. This is 
different from the synchronization method proposed by 
Pecora and Carroll [1] where only a finite number of 
possible couplings exists, that is given by the number of 
(stable) subsystems of the dynamical system. 

What makes the above introduced method for con- 
structing synchronized systems interesting for applica- 


(6) 


tions is the fact that in many cases the function s(t) 
can be rather general. In particular, it may depend not 
only on the state x but also on some information signal 
i(t), i.e., s = h(x,i) or s = h(x,s,i). This feature can, 
for example, be used in a communication scheme where 
x = f(x,s) is the (chaotic) dynamical system of the trans- 
mitter, s = h(x,i) is the transmitted and received signal, 
and y = f(y,s) constitutes the receiver. In the case of 
synchronization, i.e., when y — x, the information i can 
be recovered without error from s = h(x,i) = h(y,i) if 
this equation is uniquely solvable for i. 

This encoding method has been implemented experi- 
mentally using an APD of Chua’s circuit [15]. The trans- 
mitter and the receiver are given by two copies of the 
following system: 


dV, 
Ci ie = G(Ve2 — Ver) — g(Vc1) — Gals — Ver), 
dv ; 
Ce EVs n, (7) 
dir 
Ok = Vo. 
dr c2 


The transmitted signal s = Vc; + i is given by the 
voltage Vcı at the capacitor Cı of the transmitter 
and the information signal i. The parameters are 
Cı = 10 nF, Cy = 100 nF, L = 18 mH, G = 1/1700 Q, 
a = 3.4 and g is a piecewise-linear function defined by 
g(V) = moV + (m; — mo) [IV + B,| — |V — Bpl] with 
mo = —0.409 mS, mı = —0.756 mS, and B, = 1.08 V. 
The circuit diagram and other details of this implemen- 
tation will be given elsewhere [16]. Figure 1 shows as 
an example the reconstruction of a triangular informa- 
tion signal i(t). Although the parameters of the two 
coupled circuits are not exactly the same, the quality 
of the reconstructed signal [Fig. 1(b)] is already quite 
good. In general this sensitivity of the synchronization 
on parameter differences may cause difficulties for typical 
hardware implementations. On the other hand, however, 
it is a feature that is very welcome for any private com- 
munication. This problem and the influence of additional 
noise in the transmission channel will be discussed in 
more detail elsewhere. 

Of course, the APDs of the Lorenz system introduced 
above may also be used for encoding. System (5), for 
example, can be driven by s = 10x2 + ix3. Note that 
the above given proof for the synchronization of this 
system (e — 0) holds for all driving signals s. Therefore, 
y converges to x for all information signals i (that, of 
course, have an effect like dynamic noise and thus modify 
the underlying attractor). Since the variable x3 is always 
positive, the information signal can be recovered exactly 
as ig = (s — 10y2)/y3 as soon as the transient of the 
synchronization process is over. For the second system 
(6) the transmitted signal s may be generated by the 
differential equation s = 10(x; — s) + i. To retrieve the 
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FIG. 1. Experimental implementation of a communication 
scheme using Chua’s circuit (7) with driving s = Vc; + i and 
a triangular information signal i. (a) Transmitted signal s. 
(b) Recovered signal iz. 


information in the receiver as ig = s — 10(y, — s) in this 
case the derivative s has to be computed from s. Of 
course, other functions h (including linear filters) could 
also be used to generate the transmitted signal s. The 
only restriction for the signal s is that it has to be chosen 
in a way that the transmitter and the receiver remain stable 
(and chaotic). 

The last example shows that the Pecora-Carroll syn- 
chronization method can also be used in a more sophis- 
ticated way for communication than it was done until 
now. The communication scheme presented in Refs. [2,3] 
is based on Pecora-Carroll synchronization and cascaded 
subsystems in the receiver that are necessary to recover 
the information signal. Compared to the method dis- 
cussed in this Letter the main difference consists in the 
fact that in Refs. [2,3] the information is just added to a 
chaotic carrier but not injected into the dynamical system 
constituting the transmitter. In this case, the receiver is 
forced by a sum of the chaotic signal and the information 
signal whereas the transmitter is just driven by the (pure) 
chaotic signal. Because of this (slightly) different driv- 
ing y does not converge exactly to x and the information 
signal can only be recovered including an error that van- 
ishes in the limit |i] — 0. The specific properties of this 
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inevitable error are investigated in Refs. [4,5]. If one uses 
an information signal with a small amplitude |i| in order 
to minimize this error, however, the information can be 
destroyed by noise in the transmission channel. Further- 
more, this method is not very secure, because it is pos- 
sible to fit a nonlinear model to the time series given by 
the transmitted signal s that consists mainly of the (low- 
dimensional) chaotic carrier (|i| has to be small to avoid 
errors). Using this nonlinear model the information sig- 
nal may then be extracted from s by methods similar to 
nonlinear noise reduction [17] or other techniques [18]. 
In contrast, even if we use for our communication scheme 
just a linear combination like, for example, s = x. + i, it 
is difficult to decode the information from the transmit- 
ted signal s since the dynamical system of the transmitter 
is not autonomous but driven by the information signal 
which is in general rather complicated. 

The communication method discussed above can also be 
implemented using high-dimensional systems and trans- 
mitted signals that are hyperchaotic. In this case, how- 
ever, it is in general difficult to find a region in parameter 
space where hyperchaos exists. Therefore, we propose an- 
other method for synthesizing high-dimensional systems 
in a systematic way, using standard low-dimensional sys- 
tems with a well-known dynamics, only. For the sake of 
brevity, only the procedure for designing a six-dimensional 
transmitter using two three-dimensional systems as build- 
ing blocks is described in the following. The equations of 
the transmitter, the transmitted signal, and the receiver are 
in this case given by 


x} = fi (X1, Saux(t)) 
Saux T h(x, i(t)) 
x2 = f2(x2, s(t)) 


transmitter, 


s = hy(X2, Saux) transmitted signal, 


= foly2, s(t)) 
Sad i hy '(y2, s(t)) 
4 = fi(y, Saux (t)) 


where we assume that h; and hz are invertible with respect 
to i and Saux, respectively. If both pairs of subsystems 
(x1, yi) and (x2, y2) synchronize mutually then at the re- 
ceiver the information ig can be recovered as 


ir = Ay (Yi. Sau) = Ay (yi. hy '(y2,8))- 
The generalization to a communication model with an ar- 
bitrary number of low-dimensional subsystems is straight- 
forward [11]. The two low-dynamical systems we use in 
the following numerical example are the Réssler system 
and the Lorenz system. The transmitter of our communi- 
cation model is given by 


x1 =2 + x1 (x2 a 4), 


receiver, 


x4 = — 10x4 +s r 
X2 = —x1 — X3, X5 = 28x4 — x5 — xX4X6, 
X6 = X4X5 — 2.666x6 , 


s= 10x5 + 30Saux/X6 > 


x3 = XO 2.45.x3 + Saux> 


Saux = 1 + 3x3, 
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where s is the transmitted signal. The receiver reads 
yı =2+ yidr — 4), 
ys = 28y4 — ys — yao, Y2 = TYI T Y3, 


y4 = —10y4 + s, 


Yo = yays — 2.666y6, y3 = y2 — 2.45y3 + Saux » 
Saux = (s — 10ys)y6/30, ir = (s — 10ys)y6/30 — 3y3. 


The parameters of the Rössler and the Lorenz equa- 
tions are chosen such that both systems have a chaotic 
attractor. The information signal shown in Fig. 2(a) is 
the spoken word “42” recorded with a sampling rate of 
8000 Hz and a resolution of 16 bit. Figure 2(b) shows 
the transmitted signal s. The transmitter and receiver syn- 
chronize, and the difference |i — ig| between the origi- 
nal information signal i and the reconstructed signal ip 
is given in Fig. 2(c). The transmitted signal is in this 
case a sum of two chaotic signals and hyperchaotic. 
The Lyapunov exponents of the transmitter for i = 0 
equal A, = 0.710, A2 = 0.177, Az = 0.000, Ay = —0.006, 
As = —2.59, and As = —14.37. This and other examples 
of high-dimensional systems with more than one positive 
Lyapunov exponent [11] show that it is possible to syn- 
chronize hyperchaotic systems using a scalar signal, only. 


0 20 40 60 80 100 
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FIG. 2. Numerical simulation of a communication scheme 
using a combination of Réossler and Lorenz systems. 
(a) Information signal i given by the spoken word “42.” 
(b) Transmitted signal s. (c) Difference |i — ig| between the 
original and the recovered information signal. 


In conclusion, we have presented a general approach to 
construct a large class of synchronized dynamical systems 
and discussed its application in a communication scheme 
where the information can be recovered without errors. 
Furthermore, we introduced a procedure for synthesizing 
high-dimensional synchronized systems. In this case the 
transmitted signal is hyperchaotic and therefore yields a 
more efficient encoding. 
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